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Abstract 

The nuclear structure function F%a{x) has been studied in the Bjorken 
limit for (1,1') scattering on nuclei in the region of x > 1 and was found 
to be very sensitive to the information contained in the nucleon spectral 
function in nuclei, particularly the correlations between momenta and 
energies in the region of large momenta. Calculations were done in a 
local density approximation using two different spectral functions for 
nuclear matter. Results are compared to those obtained for a spectral 
function which has been evaluated directly for the finite nucleus, 16 0, 
under consideration. For values of x around 1.5 and larger the quasi- 
particle contribution is negligible, thus stressing the sensitivity of the 
present reaction to the dynamical properties of nuclei beyond the shell 
model approach. Several approximations which are usually employed in 
studies of the EMC effect have been analyzed and their inaccuracy in 
this region is demonstrated. The results stress the fact that the nuclear 
structure function contains important information on nuclear dynamical 
correlations. Therefore further measurements of ^2,4(2;) in that region 
and for many nuclei would be most welcome. 
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1 Introduction 



Deep inelastic scattering of leptons on nuclei in the region of asymptotic free- 
dom and the ratio R(x) of the nuclear structure function compared to the 
corresponding one for the deuteron (EMC effect |TJ) has been one of the top- 
ics in the interface of nuclear and particle physics intensively studied in the 
recent past After multiple discussions, pionic effects f|,|)|-f|, binding 



effects p-O l[ and Fermi motion [Em, 12] have turned out to be important 



ingredients to describe the main characteristics of the ratio R(x): the small 
enhancement beyond unity around x ~ 0.1, the small depletion around x ~ 0.6 
and the steady increase around x ~ 0.8 and above, respectively. 

The increase of R(x) at x close to one was soon identified to be a con- 
sequence of Fermi motion |3],|6],[12| and this is one of the points where there 



seems to be consensus among scientists. Probably the lack of controversy on 
this issue prevented a systematic exploration of the region of x > 1, in spite of 
the fact that in nuclei only the variable xa = —q 2 /2MAq° is limited between 
and 1, while x = —q 2 /2Mq° varies from to Ma/M. Since F2n(x), the 
structure function of a free nucleon, is zero for x > 1, the fact that F2a{x) is 
different from zero for x > 1 must necessarily be attributed to modifications 
of the nucleon properties inside nuclei, i.e. to genuine many body effects. The 
region of x > 1 is hence a source of very interesting information on nuclear 
properties, as we shall see. 

The fact that Fermi motion was so important at x ~ 1 induced people to 



investigate the region of x > 1 with the same idea [11-14]. These works es- 
sentially employed phenomenological nuclear matter momentum distributions, 
with the main conclusion that a nonvanishing value of F2a{x) for x > l + kp/M 
(with kp the Fermi momentum) would require a tail in the momentum distri- 
bution n(k) for k > kp. The occupation number n(k) is different from zero 
for k > kp as soon as the effects of a residual NN interaction are considered 
leading to a correlated many-body system of Fermions. 

However, it is quite dangerous to reduce the effects of correlations to a dis- 
cussion of a momentum distribution only, because in a system of interacting 
Fermions the energy and momentum distributions are correlated by means of 
the spectral function, and the simultaneous consideration of both distributions 
is necessary in principle, and also in the practice of the present case, as we 
shall see. Disregarding the energy and momentum correlations leads some- 
times to quite erroneous results, like in the study of the A mesonic decay in 
nuclei where the results for the width based exclusively on the nucleon mo- 
mentum distribution in nuclei are three orders of magnitude bigger than the 
results obtained with a spectral function, or the experimental results ]15[|. The 
same warning, in the present context, was raised in ref. (16| where F2a{x) was 
evaluated for 3 He. Further work in this direction was done in refs. |T7|-[l9l. 

In the present work we have evaluated F 2 a{x) for several nuclei by using 
spectral functions of infinite nuclear matter and the local density approxima- 
tion. This approximation is good when dealing with volume processes like the 
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present one. In order to quantify uncertainties from the many body approach 
used, we have performed the calculations using two spectral functions evalu- 
ated with rather different methods. On the other hand, since large momentum 
components are necessarily involved in the process, relativistic corrections are 
bound to be relevant [0,|2l] and we have worked within a relativistic approach. 

We evaluate F 2 a{x) for values of 1 < x < 1.5. At x = 1.5 the value of the 
structure function decreases by about three orders of magnitude with respect 
to x = 1, but is still well within measurable range. We also show results using 
the momentum distributions alone, within several approximations, and find 
appreciable differences (of two or more orders of magnitude) with respect to 
the accurate results using the spectral functions. These results show that the 
study of nuclear structure functions at values of x > 1 is a very interesting tool 
to learn about nuclear dynamical correlations beyond the nuclear properties 
described in a shell model approach. 

In this paper we not only discuss results obtained within the local density 
approximation but also consider a spectral function calculated directly for the 
nucleus 16 O. After this introduction we discuss in section 2 the calculation of 
the nuclear structure function and its relation to the spectral function. The 
various approaches for the spectral functions are presented in section 3. The 
results of our studies are presented and discussed in section 4. The last section 
summarizes the main conclusions. 



2 The nuclear structure function 

Deep inelastic electron (or muon) scattering on an unpolarized nucleon can 
be described in terms of two structure functions, Wi(x, Q 2 ) , W 2 (x, Q 2 ), where 
the Bjorken variable x is given by 

x = = (1) 

2pq 2pq 

with q the momentum of the virtual photon and p the momentum of the 
nucleon. In the Bjorken limit, q° — > oo, Q 2 — > oo and x fixed, it is common to 
define the structure functions F\ and F 2 which depend only on the variable x, 
up to some smooth logarithmic dependence on Q 2 from QCD corrections. In 
this limit one has 

p ^W 2 (x,Q 2 ) = F 2 (x) 

MW^x.Q 2 ) =Fx{x) 
and F 2 (x) , Fi(x) are related by the Callan-Gross relation 

2xF 1 (x) = F 2 (x) (3) 

Using these structure functions W\ and W 2 the hadronic tensor for the 
absorption of the virtual photon can be written: 
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W ^ = + O-) + PV^ 

p" 1 = jf - P 4q» (4) 

It is practical to work in a frame where q is parallel to the z direction. 
Adopting this frame and inspecting the transversal W lxx component in the 
Bjorken limit, one finds that the term proportional to W 2 in eq. (4) vanishes 
and W' xx is related to W\ for nucleons in nuclei with the same coefficient 
in front as for one nucleon in the vacuum, independent on its momentum or 
energy. This allows one to write the structure function Fia derived from lepton 
scattering from a nucleus with baryon number A in a nonrelativistic formalism 
using eq. (2). 

Fia(x a ) f d 3 p «* F lN (x N ) 

where Sh(oo, p) denotes the hole spectral function, i.e. the probability of finding 
a nucleon with energy u and momentum p in the nucleus, the integration limit 
\x the chemical potential or Fermi energy and 

—q 2 —q 2 —q 2 
2p A q 2M A q" 2pq 
Instead of x A one normally uses the variable x, 

-q 2 M A 

X ~ 2Mq° ~ M XA U) 

so we will write Fia and F 2A as functions of x from now on. 

By means of eq. (3) for nuclear targets we can calculate F 2A which is the 
structure function used in studies of the EMC effect. We then write 

F 2A (x) = I 4^ f S h {u,p)— F 2N (x N )0(x N )6(l-x N ) (8) 
J [2Tiy i J-00 x N 

Instead of using the spectral function Sh calculated directly for the nucleus 
under consideration, it is common practice to employ a local density approx- 
imation and represent this spectral function in terms of a spectral function 
Sh(uj,p; p) evaluated for infinite nuclear matter at various densities p which is 
normalized by 

/cftp 
JM S J-~ dWSh{W ' P;P) = P ' W 
with a factor 4 on the left side of this equation to account for the spin-isospin 
degeneracy of symmetric nuclear matter. Assuming a density profile p(r) for 
the finite nucleus to be studied, one can determine the local density approxi- 
mation for the spectral function of this nucleus by 
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S h (u,p) = A / d 3 rS h (u,p;p(r)) 



(10) 



which ensures that 



4 / d 3 r 



d 3 p 
(2^)3 



do; (u,p;p(r)) 



d 3 p 
(2^)3 



do; Sh{oJ,p) = A . (11) 



For nuclear matter the spectral function can be evaluated in terms of the 
nucleon selfenergy T,(u,p) by 



ImYi(tju,p) 



(12) 



7T [u - e(p) - ReZ(uj,p)} 2 + [/mS(u,p)] 2 

where we have dropped the variable p identifying the density dependence of 
the self-energy and spectral function. In eq. (12), e(p) is used to represent the 
nucleon kinetic energy. 

Relativistic corrections accounting for the kinematics of the nucleon have 
been included in deep inelastic scattering. Prescriptions based on the normal- 
ization of the relativistic current operator lead to corrections of the static, or 
shell model, structure function of the nucleus |2"2"| . Further corrections have 
been considered in |20|.|21~1 



23 a different relativistic treat- 



at x ~ 1. In ref. 

ment is developed which allows to write all quantities in terms of the nucleon 
propagators. Only the region < x < 1, which was measured by the EMC 
collaboration, is studied there. As in ref. 0,^,0], pionic corrections are shown 
to be relevant in the region of x < 0.6, but they play no role in the region 
x > 1 which we study here. 



Employing the treatment described in ref. [23], which uses a relativistic 



spectral function from the beginning, one can avoid introducing any flux fac- 
tors as in ref. 



22 



to account for relativistic corrections in non-relativistic 
nuclear wave functions. 

Since the relativistic corrections are important here we take advantage to 



discuss briefly and complement the details of ref. [^3 



In fig. la we show the Feynman diagram that symbolizes the deep inelas- 
tic process on a nucleon. The final hadronic state X will contain at least one 
baryon and will have baryonic number one. In a nucleus the nucleon N will 
have a certain momentum and energy distribution given by the spectral func- 
tion. The most practical way to take this into account and have a covariant 
formulation of the nuclear problem is to fold the amplitude in fig. la and 
convert it into a many body diagram for the selfenergy of an electron in the 
nuclear medium, fig. lb. Here the nucleon N in fig. la gets converted into a 
hole line and with the baryon existing in X it completes a fermionic loop. In 
section 3 of ref. [23] one evaluates this selfenergy in infinite nuclear matter and, 
by means of the local density approximation, the (e, e') cross section, is related 
to the imaginary part of the electron selfenergy. The imaginary part of the e 
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selfenergy is evaluated using Cutkosky rules and this means that the interme- 
diate states e' and X are placed on shell in the integrations over the momenta 
of these states. The formalism is originally covariant in the sense that every- 
thing is written in terms of propagators of the particles and we can write the 
nucleon propagator in a covariant relativistic way. However, Cutkosky rules 
select the imaginary part of the propagator of nucleon N for the occupied 
states and when doing that the apparent covariant structure might not show 
up clearly. 

In our formalism we start from a free nucleon propagator which we split 



(13) 



into its positive and negative energy parts p4 

rf + M = M f T, r U r (p)U r (p) Er Vr(-p)Vr(-P) \ 

p 2 — M 2 + ie ~ E{p) \p° — E(p) + ie p° + E(p) - ie J 

where M, E(p) are the nucleon mass and the relativistic nucleon energy (p 2 + 
M 2 ) 1 / 2 and u r (p) , v r (p) are the ordinary spinors which we take normalized as 
u r (p)u r (p) = 1. We recall that u r (f ) are functions of three momentum and 
they are the only spinors which appear in our framework. 

In order to account for binding and momentum distribution of the occupied 
nucleons we need the nucleon propagator in the nucleon medium. 

Note that even if a nucleon is off shell, p° ^ E(p), in the propagator of 
eq. (13) and we have p 1 + M — > p ^f — fry + M, the positive energy part has 
the Dirac structure 2M J2r u r (p)u r {p) = E(p)j° — + M, corresponding to 
on shell nucleons of momentum p. 



Following a standard relativistic notation |25|] the nucleon propagator in a 
spin saturated system would be 

G (P°.P) = ^-M-E.-E*/ (14 > 

which includes a scalar and vector terms in the nucleon selfenergy (the inclu- 
sion of a term of the type 7p does not change the arguments and conclusions 
which follow). The extraction of hole and particle spectral functions requires 
the evaluation of ReH s,v and particularly JmS s '' !; , which is a non trivial task. 

We respect the structure of eq. (14) but follow a different path in order 
to single out the imaginary part of the positive energy piece of the nucleon 
propagator. We start from the realization that for this latter purpose, in 
a perturbative expansion of the propagator of eq. (14) in terms of the free 
propagator of eq. (13), the terms of positive energy will be singular and 
dominate over those of negative energy. This allows us to write the desired 
part of the propagator as 

M ^ u r (p)u r (p) M u s (p)u s (p) 

^ r.0 _ IP. ({ft W ' V > TP. Oft ^ ^0 _ TP.Oft + "' 



E(p)^p°-E(p) ^ E(p) V P° - E(p) 
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M x ^ u r (p)u r (f) 



(15) 



M 



E(p) 



This expansion is rather useful because both S s and E^7°( and 7*p) are 
diagonal in the base of the u r (p) spinors, which converts eq. (15) in an ordinary 
geometric series, not a matricial series, which can be summed trivially as shown 
in the last step of eq. (15). 

It might look surprising that one obtains a Dirac structure u r (p)u(p) in 
G(p°,p) as for free nucleons, even when the renormalized nucleons will be off 
shell. This is less striking if one recalls that also in eq. (13) the positive 
energy part (corresponding to G) has the same structure u r (p)u r (p) even if 
the nucleon is off shell. In any case the structure might not match the one 
coming from eq. (14) and one has lost the covariance shown by eq. (14). 
The reason for this loss of covariance is that one loses terms with admixture 
of the positive and negative parts of the nucleon propagator of eq. (13) in 
the perturbative exapansion. This is, in the sum of eq. (15) one is summing 
terms of the type of fig. 2a,b,c,d, where in the intermediate fermion lines one 
only has the part of positive energy of the propagator. One is missing terms 
of the type of fig. 2e (where the line pointing down stands for the negative 
energy part of the free nucleon propagator), which would naturally appear in 
a covariant expansion of eq. (14). 

On the other hand while all the terms of fig. 2c, d etc. are summed up 
automatically in eq. (15) in terms of a selfenergy given exclusively by the 
term in fig. 2b, the second order terms of fig. 2e is not accounted for. We 
argued that this latter term (which has an intermediate propagator of order 
1/2M, and is of p 2 type) should be small with respect to the diagrams contained 
in figs. 2b, c, d, ... But in any case it can be included as a nucleon selfenergy 
part in the sum of eq. (15) and then diagrams e,f, etc. would be automatically 
included. This means that even if the covariance is lost in eq. (15) one can 
still regain all the terms in the series by including these mixed terms in the 
selfenergy £ appearing in eq. (15). Of course, this selfenergy is now different 
to the one appearing in eq. (14). These mixed terms are also diagonal in u r (p) 
and do not change the structure of eq. (15). This is the philosophy which we 
follow, only that the diagrams of fig. 2e, 2f are not evaluated, although they 
are implicitly accounted for as we pass to discuss. The reason is that these 
diagrams only contribute to the real part of E, not to the imaginary part, and 
in our scheme, which evaluates accurately JmS, there are pieces missing in the 
real part of S which are added phenomenologically in order to ensure the exact 
experimental binding energy of each nucleus f23]| . The particular structure of 
eq. (15) allows one to write 
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+ f^sM^L. (16) 



with the relationships 

M 



ImE(p ,p) 



Sh(p°,p) 



1 E(p) 



Sh(p°,p) 



n \p° - E(p) - -^RenP^p)} 2 + [^) /mS (P°^)] 2 

for p° < fi 



-imX(p°,p) 



.1 E(pY 

* [p° - E(p) - JLjReX(p°,p)] 2 + [^/mS(p°,p)]2 

for p° > /j, 



(17) 



k F , p (p) = [3tt 2 Pp (p)] 1/3 k F>n (f) = [37r 2 P n(r-)] 1/3 



where for simplicity S is now mSu which is independent of spin. 

By means of this new nucleon propagator the modifications introduced by 



our relativistic formalism, described in detail in ref. [23|, are rather intuitive, 
easy to employ and can be summarized as: 

i) The normalization of the spectral function which ensures the proper 
normalizations of the charge (or baryonic charge) of a nucleus is exactly 
the same as in eq.(9). However, the spectral function of eq.(12) is now 
replaced by eq. (17) 

ii) On the other hand the structure function F 2 a(x) of eq.(8) is replaced by 



F 2A (x)=A[ d 3 r j ttt^t^t f du S h (cu,p)—F 2N (x N )e(x N )6(l-x N ) 

J J (2TT) 6 E[p) J~oo Xn 

(19) 

where the relativistic factor -M=r plays the role of a Lorentz contraction factor, 
appearing in the probability per unit time of electron collision with the nucleon, 
and which remains in the formula of the nuclear cross section because one 
divides the sum of all probabilities by a unique electron flux, the one of the 
electron with respect to the CM of the nucleus. 

The questions of normalization and conservation of baryonic number, which 
have been the subject of much attention [Q, are discussed in detail in |25fl . 

There are other terms which would be included in a covariant formalism 
of the (e, e') reaction and do not appear in our formalism. These are terms 
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which have a negative energy state coupled to the hadronization vertex, as 
shown in fig. 3. Once again such terms are reduced by the large energy 
denominator of the negative energy state. Only in cases when one uses an 
operator which magnifies the NN coupling with respect to the NN, as in 



the case of the axial charge, such terms can be relevant |26|, although more 



accurate nonperturbative calculations find smaller results P7fl . With the use 
of scalar potentials S s smaller than the typical ones in the Walecka model 
when one imposes constraints from information in the negative energy sector, 
as done in [p8|| , terms like in fig. 3 would be of the order of 10-15 % if one 
has an operator like in the axial charge, 7°75, but much smaller than this if 
electromagnetic current operators are used ||27| . 

Furthermore, in as much as one assumes that the structure functions for the 
positive or negative energy states are the same and one uses a nucleon spectral 
function which conserves the baryonic number, one would be including the 
strength of these pieces into the scheme which we follow, up to small differences 
coming from different medium corrections to the positive and negative energy 
states. Estimates based on the findings of [T27J] would put the difference between 
this covariant scheme and ours at the level of 1-2% in the EMC region and 
probably a few percent in the x > 1 region that we explore. 

In practical terms our scheme respects special relativity in the positive 
energy sector and amounts to using free u r (p) spinor in the evaluation of the 
matrix elements of the 7*iV — > X process while keeping the proper energy and 
momentum balance in the 5 function of conservation of fourmomentum, with 
the u,p distribution of the occupied nucleon given by the spectral function 
and the energy and momentum of the final states in X being those of their 
asymptotic states. This is in fact the most standard method in the study 
of many nuclear processes involving scattering or decay. The aproximations 
which we have done here, sacrificing covariance in a controled way, lead us to 
this calculational scheme where everything is defined. Covariant formalisms 
like those used in H29], |3~0l generate some off shell dependence in the hadronic 
tensor, which are accounted for in terms of new structure functions for which 
there is no empirical information, so several different assumptions are made 
in |?I||3l| which produce moderate changes in the EMC results. 

Our relativistic corrections thus stem from the consideration of special 
relativity in the positive energy sector, although, as we discussed, it accounts 
in an approximate way for the contributions involving negative energy states, 
which are small anyway. 

We will show results both with the relativistic and non-relativistic formal- 
ism. The relativistic corrections are found to be relevant in the region of x > 1, 
particularly at large values of x. 

It is easy to see qualitatively why eq.(8) or eq.(19) lead to a non- vanishing 
structure function for x > 1. The nucleon structure function appears with 
argument xjy in this equation and in the Bjorken limit one has 
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%n = — x = x (20) 

X uj — p z 

For certain combinations of uj and p z one can obtain values of xn < 1 even 
if x > 1. Since both uj and p z appear in eq.(20) it is very important to take into 
account the correlations between uj and p provided by the spectral function 
Sh(uj,p), and one sees that approximations which neglect these correlations 
are bound to provide unrealistic results. 

3 The nucleon spectral function in nuclear mat- 
ter and finite nuclei 

We have used three different approaches to evaluate the spectral function. 
The first one is a semiphenomenological one relating the spectral function 
of nuclear matter to the experimental cross section for NN scattering. The 
second approach is microscopic in the sense that the spectral function for 
nuclear matter is derived from a many-body calculation employing a realistic 
One-Boson-Exchange model for the NN interaction. In the third model we 
avoid the local density approximation and evaluate the spectral function Sh 
directly for finite nuclei. We briefly describe these models below. 



3.1 Semiphenomenological approach 



This model is described in detail in ref. |JT|]. It evaluates /mE(w,p) from a sec- 
ond order diagram and uses the fact that ladder diagrams evaluated from the 
iViV potential lead to the NN t matrix. Pauli blocking corrections are taken 
into account in the explicit diagram evaluated and |t| 2 , which appears in the 
evaluation of JmS, is written in terms of the experimental NN cross section. 
Polarization effects from the RPA iteration of ph and Ah excitations are also 
taken into account. The real part of the selfenergy is obtained via a disper- 
sion relation and the Fock term from pion exchange is also included. Hartree 
terms, which require the explicit knowledge of a potential, are missing in the 
approach but these are terms independent of energy and momentum. The 
nucleon properties evaluated in ]3l] as a function of uj — fi compare favourably 



with those of more microscopic evaluations [|32],[33|]. In order to complete the 
model and obtain absolute values for i?eS, another phenomenological piece is 
added here. This Hartree contribution is assumed to be proportional to p and 
its value is adjusted, in order to fit the empirical value of the binding energy 
per nucleon in each particular nucleus. For this purpose we recall the sum rule 



for the binding energy per nucleon [34 



\e A \ = - l -( <E -M> ^ < T > ) (21) 



and we evaluate < T > and < E > as 
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By means of this, one takes also into account empirically contributions 
from terms like in fig. 2e, which are in principle small, and even if they are of 
p 2 type would not differ appreciably from a pp e // form, with p e ff an effective 
average nuclear density. 

The evaluation of ImY,(uj,p) was done nonrelativistically in [31]. For con- 



sistency with the relativistic formalism used here we should have kept the 
factors ^ in the nucleon propagators evaluating ImE(uj,p) in ref. [[Hj]. How- 



ever, the range of momenta in the loop integrals in /mE(w,p) is quite limited 
and they would modify the values of JmE(w, p) by less than 10%. By means of 
the explicit calculations carried out here we have observed that an increase of 
10 % in ImT>(u, k) leads to increases of F 2 a(x) of the order 2 % at < x < 0.6 
and always below the 10 % level for large x, hence we have continued to use 
the same ImH(uj,p) as obtained in ref. fl3T |. 



3.2 Microscopic approach in nuclear matter 

The spectral function of nuclear matter which has been used in this sec- 



ond approach has been evaluated using the techniques described in ref. [[35 
The starting point of this many-body calculation is a Brueckner-Hartree- 
Fock (BHF) calculation of nuclear matter considering the realistic One-Boson- 



Exchange (OBE) potential B as defined in [|36| for the NN interaction. The 
G-matrix resulting from this BHF calculation as well as the BHF single-particle 
spectrum €bhf{p) are used to define the nucleon self-energy including all terms 
up to second order in G. The single-particle Green's function g(p, uj) is derived 
from a solution of the Dyson equation 

g(p, u) = g( BHF \p, lj) + g (BHF \p, uj) p 2 \ P , co)] g(p, u) (23) 

Here S^ 2 ^ is the contributions to the self-energy in second order. The 
single-particle Green's function in the BHF approximation is given by 

g(BHF) M= Q(k F -p) + Q(p-k F ) 

to — dp — it] uj — e p + ir] 

where kp denotes the Fermi momentum of nuclear matter at the density under 
consideration. The term with S^ 2 ^ in eq. (23) contains a contribution with 
intermediate particle-particle states, which is taken into account already in the 
BHF aproximation. This doublecounting is removed as described in |35| . The 
spectral function Sh(u,p) can then be calculated from the imaginary part of 
the single-particle Green's function by 
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S h (uj,p) = -Img(p,tj) (25) 

7T 

This calculation yields a Fermi energy \i depending on the density of nuclear 
matter. The energy variable u is defined with respect to this Fermi energy. 
In the local density approximation for the spectral function discussed above, 
the empirical Fermi energy of the finite nucleus has been chosen to be the 
reference point for the energy variable u. 

3.3 Microscopic approach for finite nuclei 

The spectral function can be calculated directly for finite nuclei using the 
procedure described and applied to 16 O in ref. f37|. For nuclei with spheri- 
cal symmetry the self-energy is evaluated in a partial wave basis, Ey(p,p'), 
assuming that orbital angular momentum I and total angular momentum j 
are conserved quantum numbers. As discussed above, the total self-energy is 
decomposed in a BHF part and terms of second order in the Brueckner G- 
matrix. The corresponding single-particle Green's function can be evaluated 
by solving a Dyson equation of the form 



9i,j(p,p'; w) = gij HF) (p,p'; w)+ / dk x J dk 2 gij HF \p, h; u) 

SyOi, h; w)l g hj (k 2 ,p'; u) (26) 



x 



The spectral function for the various partial waves is then obtained from the 
imaginary part of the Green's function gij(p, p; to) applying eq.(25). A problem 
of this partial wave expansion for the momentum distribution is related to 
the fact that non-negligible contributions are obtained at large momenta and 
energies in high partial waves. Therefore we prefer to apply an approach 
which has been introduced and discussed in ref. ||35|| . In this approximation 
one splits the spectral function, for nuclear matter as well as for finite nuclei, 
into a quasiparticle contribution describing the contribution to the spectral 
function around the quasiparticle pole and a background contribution which 
contains the information about the spectral function at energies away from the 
respective quasiparticle pole. For finite nuclei a quasiparticle pole contribution 
is only observed for those partial waves, which are occupied in the HF or 
independent particle model. Therefore the sum on partial waves in 

S® p (uj,p) = £2(2j + l)n ld 5( U - eff)\^(p)\ 2 (27) 

is restricted to Z=0 and 1 in our example of 16 0. In this equation e®J stands 
for the energy of the quasiparticle pole, $ij(p) for the corresponding single- 
particle wave function in momentum space and nij for the occupation prob- 
ability for this pole. This quasiparticle pole contribution is supplemented by 
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the background contribution calculated in a local density approximation from 
the background contribution in nuclear matter 

S h>A (uj,p) = S QP (uj,p)+4 J d 3 rS*(uj,p;p(r)) (28) 

where stands for the background contribution of the spectral function 
calculated in nuclear matter at the local density p(r). Care is taken that the 
whole spectral function is normalized such that 

!§?I-^ SM=A (29) 

with A the total number of nucleons. 



3.4 Approximations to be avoided 

If no reliable model for the spectral function for nucleons in nuclear mat- 
ter is available, one may be tempted to use certain approximations. One of 



such approximations, which has frequently been used fll2|-|i"4"[ is to ignore the 
special correlations between momentum and energies of nucleons provided by 
the spectral function and simply use the energy-integrated spectral function, 
which is the momentum distribution. We are going to discuss three differ- 
ent approximations and try to explore their reliability by comparing with the 
results obtained in the more sophisticated models for the spectral function 
discussed above. Although discussions around different approximations to the 
nuclear wave functions, and other different approximations, have been com- 
mon in the past [|| f|,[7],|38l, the comparison of different approximations to the 
results obtained using spectral functions has not been exploited, particularly 
in the region of x > 1 which we study here. 

a) The uncorrelated Fermi sea distribution. 

This is the simplest approximation, which is usually very accurate, except of 
course in processes which test the momentum distribution at momenta which 
are large compared to the Fermi momentum, as is the case in the present 
problem. In this approximation the spectral function within the local density 
approximation is assumed to have the form 

S% FS (u,p; p) = n FS tf)6(u ~ E(P) ~ S) (30) 
with an occupation probability of 

n F s(p)-{ ifb1>pF ( r) ( 31 ) 

with a local Fermi momentum p^(r) which is related to the local density p(r) 
by 
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p,(r) - (^) 1/3 (32) 
and an expression for the local single-particle potential 

£ = S(r) = V TF (r) + Dp{r) (33) 

where Vtf{ t ) is the Thomas Fermi potential, —p 2 F (r)/2M and .D a phenomeno- 
logical constant fitted to reproduce the binding energy per nucleon in the 
nucleus, as done in section 3.1. 

b) Use of the momentum distribution of the correlated Fermi sea. 

Since large momentum components are needed to generate F 2 a(x) at x > 1, 
one is tempted to use the realistic momentum distribution of the nucleus as a 
way to improve on this approximation. This means that we assume an expres- 
sion for the spectral function Sff D like in eq. (|30|) but replace the momentum 
distribution of the free Fermi gas nps by the momentum distribution of the 
interacting Fermi gas 

n I (p) = / S h (u,p)du; (34) 

J -co 

where we have used the spectral function of nuclear matter discussed above 
to calculate nj. Note, however, that the energy- momentum relation is still 
determined by the 5-function in (30), with £ defined with the same prescription 
as in the subsection above, eq. (33). 

c) Use of the correlated momentum distribution and the corresponding mean 
value for the energy. 

Finally we want to consider an approximation in which we assume again 
a definite relation between momentum and energy of a nucleon in the hole 
spectral function 

Sff ED (^,P;p) = n z (p)6{u>- < u(p) >) (35) 

but determine the momentum distribution nj (see eq.(34)) as well as the mean 
value of the energy for a given momentum 

<u{p)>= — — 36 

JZaoShfapjdU) 

from the complete spectral function for nuclear matter discussed above. 



4 Results and discussion 

In a first step we want to compare the two approaches to determine the spectral 
function for nuclear matter, which we describe in sections 3.1 and 3.2 and which 
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we are going to employ for the calculation of the structure function. For that 
purpose we present in Fig. 4 the momentum distribution nj(p) (see eq.(34)) 
calculated at the empirical saturation density, p = p , of nuclear matter. 

The momentum distributions obtained by these two very different methods 
are very similar. At small values of p the microscopic approach of |35[] provides 
an occupation number of the order of 3 % bigger than the semiphenomeno- 



logical one |31| . For momenta above the Fermi momentum the distributions 
are also similar although for momenta around two times the Fermi momen- 
tum (~ 550 MeV/c) the differences become more appreciable. The semiphe- 
nomenological approach provides a little less strength below the Fermi mo- 
mentum, which is then redistributed to larger momenta where ni(p) is larger 
than in the microscopic approach. The precise values of ni(p) calculated in 
a microscopic many-body theory, depend on the model of the NN interaction 
which is considered and the method which is used to determine the effects 
of correlations. For instance in the self-consistent Green's function approach 



of ref . p3| , using the Reid soft-core potential, the occupation number for mo- 
menta below kp is around 0.85, smaller than in both the approaches considered 
here. This demonstrates that the semiphenomenological approach provides 
a result which is in agreement with the microscopic calculations within the 
uncertainties of the microscopic approach caused by the approximation in the 
many-body theory as well as NN interaction. The differences found in rij(p) in 
the two approaches discussed have little repercussion in the values of F 2 a(x), 
which come very close to each other in the two approaches, as we shall see 
below. 

As a second quantity characterizing the bulk properties of the spectral 
functions calculated by these two methods, we show in Fig. 5 the mean value 
for the energy as a function of p calculated according to eq.(36) at the nuclear 
density po- Fig- 5 shows in a qualitative way that there is an important corre- 
lation between the momenta and the mean value of the energy for the bound 
nucleons. The absolute value of this mean energy | < uj(p) — M > | decreases 
as a function of momentum with increasing momenta for momenta below the 
Fermi momentum. This momentum dependence is mainly due to the mo- 
mentum of the quasiparticle peak, which is approaching the Fermi energy for 
p —>■ pp. There is no quasiparticle contribution to the hole-spectral function 
Sh for momenta larger than the Fermi momentum pp. Therefore at these mo- 
menta, the mean value is determined only from the background contribution. 
The coupling to 2 hole-1 particle and more complicated configurations with 
total momenta p, described by these background terms yields a mean value of 
< uj — M >, which decreases with increasing momentum. From this figure it 
is evident, however, that the energy- momentum relation obtained from a re- 
alistic spectral function is quite different from the simple relation used in eqs. 
(30), (33), which provides a dispersion relation which is always an increasing 
function with increasing momentum. 

The values obtained for | < uj{p) > | in the two approaches discussed in 
sections 3.1, 3.2 are very similar, with differences of the order of 10 % at most. 
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This is another indication that the basic features of the spectral function are 
not very sensitive to the method used in the evaluation and that also the 
semiphenomeno logical approach yields quite a reliable result. 

In Fig. 6 we show the results for F 2 a(x) calculated with the three different 
spectral functions introduced in sections 3.1 to 3.3, for the case of 16 0. The 
density distributions p(r) for 16 and the other nuclei, which are required to 
apply the local density approximation, are taken from refs. |39|,[|0]]. Since the 
microscopic nuclear matter and finite nuclei approaches are nonrelativistic, 
we have also taken the nonrelativistic version of the semiphenomenological 
approach, for comparison, omitting all the factors in eqs. (17), (19). The 
experimental values for F 2 n{xn) are taken from ref. pEI] . The results obtained 
with the two spectral functions of nuclear matter (solid line and dashed line) 
are rather similar. At x ~ 1 the microscopic spectral function provides results 
about 20 % higher than the semiphenomenological one. At values of x ~ 1.22 
the two approaches coincide and for x ~ 1.5, where the structure function has 
decreased three orders of magnitude with respect to the value at x — 1, the 
semiphenomenological approach provides values of F 2 a about 40 % larger than 
the microscopic one. This reflects the fact that the former model provides a 
larger probability for the momentum distribution at high momenta than the 
latter one, as seen in Fig. 4. 

The results for the structure function obtained with the spectral function 
of eq.(28) evaluated directly for the finite nucleus are represented by the dot- 
dashed line in Fig. 6. They should be compared with those displayed by 
the solid line since the background contribution to eq. (|28|) is obtained from 
the same nuclear matter result. These two results can hardly be separated 
on the logarithmic scale of the figure. We observe that at x ~ 1 the results 
with the spectral function of the finite nucleus are about 8 % bigger than with 
the nuclear matter approach. The differences become smaller as x increases 
and for values of x ~ 1.5 the two approaches give the same results. This 
latter fact is telling us that at large values of x one is getting practically all 
contributions from the background part of the spectral function and none from 
the quasiparticle part. The comparison of these two curves also tells us that 
the use of the nuclear matter spectral functions, together with the local density 
approximation, is an excellent tool to evaluate F 2 a{x). If one compares the 
results at values of x studied in the EMC effect, < x < 0.6, the differences 
among the three calculations are of the order of 3 %. 

In Fig. 7 we show results obtained with the semiphenomenological ap- 
proach using the relativistic and nonrelativistic formalisms. The trend of the 
results is similar, however, the relativistic corrections induce a reduction of 25 
% around x = 1 and roughly reduce the structure function F 2 a to one half of 
the non relativistic results at x ~ 1.5. The relativistic effects are significant 
in the sense that they are bigger than the differences found between various 
nonrelativistic approaches, which reflect the uncertainties in the treatment of 
correlations. 
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Results obtained for F 2 a{x) using the different approximations discussed 
in section 3.4 are displayed in Fig. 8. The first one, which originates from 
the assumption of an uncorrelated Fermi sea, eq. (30), is represented by the 
dot-dashed line. We can see that at x ~ 1 it already provides a structure 
function of around a factor two smaller than the one obtained with the proper 
spectral function (short dashed-line). However, as one moves to higher x, the 
discrepancies get bigger and at values of x ~ 1.2 the uncorrelated Fermi sea 
gives already values for the structure function which are about two orders of 
magnitude smaller than the correct ones. It is clear that one is exploring the 
region of large momenta, above the Fermi momentum, which are not accounted 
for by the uncorrelated Fermi sea. 

Another approximation corresponds to using the realistic momentum dis- 
tribution nj{p) of eq. (34) and associating an energy to each p given by its 
kinetic energy plus a potential, eq. (33). The results (solid line) are outra- 
geously wrong. This demonstrates that the naive use of a momentum distri- 
bution, although calculated in a realistic way, may lead to results which are 
worse than those obtained for an uncorrelated system, if one does not treat the 
energy-momentum correlation properly. As we have discussed already in Fig. 
5, the mean value of the energy < u(p) > decreases above the Fermi momen- 
tum with increasing momenta. On the other hand, the energy associated to 
p in eq. (30) grows like the kinetic energy as \p\ increases. The discrepancies 
with the correct results are about a factor three at x ~ 1 and three orders of 
magnitude at x ~ 1.5, providing a gross overestimate of the results for F 2 a(x). 
The same gross overestimate found here for this approximation was also found 
in ref. |15] in connection with the mesonic Lambda decay in nuclei. 

In view of the deficiencies of the previous approximations and the reasons 
for it, one might think that the results should be improved by replacing the 
kinetic plus potential energy, eq. (30), by the mean value of < u(jp) > cal- 
culated from the spectral function (see eq.(36). This is indeed the case (see 
curve with long dashes in Fig. 8), although the discrepancies with the exact 
results are still large enough to discourage this approximation too. We can 
see in Fig. 8 that at values of x ~ 1 (and also in the EMC region below) the 
approximation turns out to be quite good. However, for values of x ~ 1.3 and 
above the discrepancies with the correct results are already as big as one order 
of magnitude or more. 

The results discussed here stress the importance of using the spectral func- 
tion to evaluate F 2 a{x) since all the information contained in it, correlating 
energies and momenta, is very important, particularly at large x. We showed 
that some schemes which use only a partial information from the spectral 
function lead to rather inaccurate results and should thus be avoided. 

Finally in Fig. 9 we show results obtained for different nuclei. They are 
calculated at Q 2 = 5 GeV 2 . We can see that F 2 a(x)/A is very similar for 
the different nuclei. We have taken nuclei with N = Z or close by, to be 
able to use a unique Fermi sea for protons and neutrons as done in symmetric 
nuclear matter. For heavier nuclei with N ^ Z the results obtained here could 
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be easily extended by dealing with two different Fermi seas, but this would 
require the extension of the work of ref. |JT| to non symmetric nuclear matter. 
We do not expect, however, any special effects apart from those exposed here. 

When evaluating absolute values of nuclear structure function, not ratios 
to the nucleon or the deuteron, it is very important to take into account the 
Q 2 dependence of the structure function. This is particularly true for values 
of Q 2 ~ 1 — 10 GeV 2 , but even at Q 2 ~ 100 GeV 2 and above, where there is 
approximate Bjorken scaling, the Q 2 dependence is weak but one still has to 
consider it if one wants to make accurate predictions. 

For the Q 2 dependence of the nucleon structure function we have taken the 



parametrizations given in ref. W2 



It is interesting to compare our results with the scarce experimental data 
available. Our results refer exclusively to the deep inelastic contribution to 
electron nucleus scattering. The quasielastic contribution (where only one 
nucleon is knocked out in the first step eN — > e'N) is not taken into account 
in our formalism. At low values of Q 2 and x > 1, the quasielastic contribution 
is dominant |]TU[ and one has to go to values of Q 2 > 20 GeV 2 to have a 



dominance of the deep inelastic contribution |I7|. For this reason we compare 



our results with measurements done at Q 2 = 61, 85 and 150 GeV 2 in ref. [[13 



which improve the preliminary results reported in ref. [II] where much larger 
values were obtained. 

The results can be seen in fig. 10. The three theoretical curves correspond 
to each one of the values of Q 2 and the results decrease as a function of 
increasing Q 2 . The agreement with the data is qualitative. The slope as a 
function of x seems well reproduce but the theoretical results are in average 
40 % higher than experiment up to x = 1.05. At x = 1.15 and 1.3 there are 
only upper bounds which are compatible with our predictions. 

Experimental results at Q 2 < 5 GeV 2 have however, a large contamination 
of quasielastic contribution JT7HT9f. This is reflected by the large dispersion 
of the results as a function of Q 2 [f|5| and in the approximate y scaling of 



these results, which is characteristic of the quasielastic collisions. Nevertheless, 
we have also evaluated the deep inelastic contribution correspondieng to the 
results in E5|] with largest values of Q 2 . We evaluate the structure function 
corresponding to the lowest curve in fig. 1 of ref. | 45| . This corresponds to 



different values of Q 2 for each value of x since the data correspond to electron 
scattering with fixed initial electron energy (E e = 3.595 GeV), fixed scattering 
angle, 9 = 39° and variable final electron energy. We show the results in fig. 11. 
The values of Q 2 increase with increasing x. Atx=l,<5 2 = 3.11 GeV 2 and 
at x = 1.25 , Q 2 = 3.42. We can see that our results lie below the experimental 
data, particularly at large values of x. However, one can observe a tendency 
to be in agreement with the data at values of x < 0.8 if one extrapolates the 
data smoothly. In fact in fig. 1 of ref. [I5[] we see a confluence of the data for 
different values of Q 2 in the region of x = 0.4 — 0.6 with values which agree 
with our results of fig. 11. This would be in agreement with the conclusions 
reached in E5j where the large dispersion of the results as a function of Q 2 
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for large values of x indicates the dominance of the quasielastic contribution, 
while the tendency to stabilize the results at values of x < 0.8 indicate that 
this region of x is dominated by the deep inelastic contribution. In such case 
our results should be comparable to the data and this is indeed the case. 

5 Conclusions 

We have evaluated the nuclear structure function F 2 a(x) at values of x bigger 
than unity, especially in the range 1 < x < 1.5 where the values obtained are 
well within measurable range. For this purpose we have used sophisticated 
nuclear spectral functions which account for nuclear correlations and relativis- 
ts effects. The structure function decreases three orders of magnitude from 
x — 1 to x — 1.5. 

The strength of F 2 a{x) in that range of x is tied to the components of the 
nuclear wave function with nucleons of large momenta. These components are 
due to the two-nucleon correlations originating from realistic NN interactions. 
The momentum distribution of nucleons in the nuclear many-body system, 
however, is strongly correlated with the energy distribution of these nucleons. 
These are dynamical effects which go beyond the shell model picture of the 
nucleus and which are taken into account in terms of the nucleon spectral 
functions. The results for F 2 a(x) are very sensitive to the correlations between 
uj and p, to the extreme that usual approximations made in calculations of the 
EMC effect fail badly in the region of x > 1. In particular, at x ~ 1.5 the 
quasiparticle bound states (the occupied states of the shell model), which are 
only partly occupied in an interacting nucleus, give a negligible contribution 
to F 2 a{x) and all the strength comes from the background part of the spectral 
function. 

We have discussed in detail the results obtained with several approxima- 
tions which use only rough spectral functions or partial information from re- 
alistic ones, and which are often used. We showed that in this region of x 
none of them can be taken as a substitute of the calculation using the whole 
information of the spectral function. 

In order to quantify the intrinsic theoretical uncertainties of the results 
we used two different models for the spectral function evaluated in infinite 
nuclear matter and F 2 a{x) for nuclei was calculated using the local density 
approximation. A version of the spectral function for finite nuclei was used 
also for 16 0. The differences between the models were small, of the order of 
10-30 % depending on the region of x. We also found that the use of the local 
density approximation was an excellent tool, providing results very close to 
those obtained by direct evaluation for the finite nucleus. 

Relativistic effects were checked and found to be important. They re- 
duce the results for F 2A (x) obtained with the nonrelativistic approximation 
by amounts ranging from 25 % at x ~ 1 to nearly a factor two at x ~ 1.5. 

On the other hand we have evaluated F 2 a(x) for different nuclei and find 
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that F 2 a{x)/A becomes very similar for N = Z nuclei from 40 Ca on. 

The experimental results at x > 1 are scarce, particularly at large values 
of Q 2 . We compared our results with available data at Q 2 = 61, 85, 150 GeV 2 
and found our results about 40% higher than experiment, although the fall 
down with x was well reproduced. At values of Q 2 significantly smaller, Q 2 < 
4 GeV 2 , we found that our results for x > 1 where much smaller than the 
experimental data, which was in agreement with theoretical and experimental 
findings that this region is dominated by quasielastic scattering. At lower 
values of x, around x = 0.4 — 0.6 our results matched the experimental data, 
in agreement with the theoretical and experimental findings that this region 
is dominated by the deep inelastic contribution. 

The present investigation and the importance of the nucleon spectral func- 
tion for the precise determination of F 2 a{x) is telling us that measurements 
of this quantity for different values of x and a wide range of nuclei would pro- 
vide important information on the components of the nuclear wave function at 
large momenta and energies and the strong correlations between momenta and 
energy. This information would be very important as a test of the many body 
theories which are employed for the determination of the spectral function and 
would unveil interesting details on nuclear correlations, which complement our 
knowledge of nuclear structure beyond the basic information contained in the 
shell model wave functions. 
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Figure Captions: 

Fig. 1: a) Feynman diagram for deep inelastic electron- nucleon scattering and, 
b) electron selfenergy diagram associated to it. 

Fig. 2: a,b,c,d) Feynman diagrams of the Dyson series in the evaluation of 
the nucleon propagator including only intermedite positive energy states, e) 
Feynman diagram of the Dyson series with a negative energy intermediate 
state, f) Higher order terms in the Dyson series originated from the selfenergy 
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term of diagrame and its iteration through positive energy intermediate states. 

Fig. 3. Feynman diagram for deep inelastic electron scattering with an in- 
terating nucleon, through the coupling of the photon to the negative energy 
components. 

Fig. 4: Momentum distributions at p = p . Solid line: microscopic model 
of [35]; dashed line: semiphenomenological model [31]. 

Fig. 5: Mean value of the energy of nucleons as a function of p, from eq. (36), 



at p = po- Solid line: microscopic model [33[]; dashed line: semiphenomeno- 



logical model [31 



Fig. 6: Results obtained for the structure function of 16 0. Solid line: mi- 
croscopic nuclear matter model |J5] ; dot-dashed line: microscopic finite nuclei 
model eq.(23); dashed line: nonrelativistic semiphenomenological model |31~| . 

Fig. 7: Results obtained for the structure function of 16 using the semiphe- 



nomenological model PT[|. Solid line: nonrelativistic formalism; dashed line: 
relativistic formalism. 

Fig. 8: Results obtained for the structure function of 16 using different 
approximations. Dot-dashed line: uncorrelated Fermi sea, eq. (25); solid line: 
momentum distribution of the correlated Fermi sea, eq. (29); long dashed 
line: momentum distribution of the correlated Fermi sea and average energy 
< uj{p) >, eq. (30); short dashed line: spectral function, eq. (13). Q 2 = 
5 GeV 2 . 

Fig. 9: Results obtained for the structure function per nucleon for different 
nuclei. Solid line: 40 Ca; short dashed line: 56 Fe; dot-long dashed line: l2 C; 
dot-short dashed line: 16 0; long dashed line: 6 Li. All results are obtained 
using the relativistic version of the spectral function, eq. (13). Q 2 = 5 GeV 2 . 

Fig. 10: Results for the structure function of l2 C at Q 2 = 61, 85 and 150 GeV 2 
(solid, short-dashed and long dashed lines respectively). The data are from 
ref. |43| crosses for 61 GeV 2 , squares for 85 GeV 2 and triangles for 150 GeV 2 . 
The data for the two largest values of x are upper bounds. 

Fig. 11: Results for the deep inelastic structure function of 56 Fe at different 
values of Q 2 around 3 GeV 2 , see text, compared with experimental inclusive 
data results of ||45|| . The experimental data for values of x > 1 are dominated 



by the quasielastic contribution while for values of x < 0.8 the deep inelastic 
contribution dominates the reaction. 
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